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PREFACE 


This report documents the work done on NASA Grant NAG8034 
to establish the control and direction of effective eigenvalue 
excursions of lightly damped, speed dependent rotor systems 
using passive control. Both second order and sixth order 
bi-axis, quas i- 1 inear, speed dependent generic models were 
investigated. In every case a single, bi-directional control 
bearing was used in a passive feedback stabilization loop to 
resist modal destabilization above the rotor critical speed. 

Assuming incomplete state measurement, sub-optimal control 
strategies were used to define the preferred location of the 
control bearing, the most effective measurement locations, and 
the best set of control gains to extend the speed range of 
stable operation. Speed dependent control gains were found by 
Powell’s method to maximize the minimum modal damping ratio 
for the speed dependent linear model. An increase of 300 
percent in stable speed operation was obtained for the 
sixth order linear system using passive control. 

Simulations were run to examine the effectiveness of the 
linear control law on nonlinear rotor models with bearing 
deadband. The maximum level of control effort (force) 
required by the control bearing to stabilize the rotor at 
speeds above the critical was determined for the models with 


bearing deadband. 
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CHAPTER I 


INTRODUCTION 

This study investigates the stabilization of 
self-excited rotor vibrations by active control. Rotor- 
dynamic systems under the influence of self-exciting 
mechanisms exhibit unstable behavior above some thresh- 
old speed. Quas i- l inear, speed dependent models of 
these systems show eigenvalue excursions into the 
unstable zone as this threshold speed is approached. 
Nonlinear models with bearing deadband exhibit unstable 
oscillations above this same threshold speed. By the 
application of control forces as can be produced by 
magnetic bearings, the eigenvalues of the linear models 
are re-positioned such that the systems are stable. The 
control is then applied to the linear and nonlinear 
models by simulation. The effectiveness of the control 
is judged by the resulting stability characteristics of 
the nonlinear models and by the control bearing force 
levels required to maintain rotor stability. 

Two rotor models are used in this study. The 
first of these is a two degree-of -freedom model resem- 
bling a modified Jeffcott model. Control bearing forces 
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are structured to counteract the self-exc iting forces, 
resulting in stable operation at any rotor spin speed. 

By the use of a speed dependent control strategy, the 
eigenvalues of the closed-loop system are maintained at 
fixed locations. The control is chosen such that these 
locations match those of the original system at its 
optimally damped speed. 

The second rotor model used in this study is a six 
degree-of -freedom model. It is shown to possess three 

natural modes, each becoming unstable at a different 
rotor spin speed. The stabilization of all three modes 
using one magnetic bearing set is investigated. The 
ability of the control bearing to stabilize each mode is 
dependent upon such factors as bearing location, feed- 
back signals, and control gains. The magnetic bearing 
location is determined by examining the uncontrolled 
rotor’s mode shapes, with placement attempting to maxi- 
mize the bearing's effect on each mode. Feedback is 
provided through combined output signals with incomplete 
state information. The combination of the output sig- 
nals and the control gains are determined by using a 
search method, which maximizes the damping of the least 
damped mode. Control strategies using constant and 
speed dependent control gains are compared based upon 
the effect that each has on the excursions of the linear 
systems’ eigenvalues towards the unstable zone. 


CHAPTER II 


BACKGROUND 

The first published work dealing with the dynamics 
of rotating shafts was by Rankine [1] in 1869. At that 
time Rankine showed that rotating systems were unstable 
above the critical speed (the speed at which synchronous 
resonance occurs). Because of this, for many years 
manufacturers designed rotating systems to operate below 
the critical speed. In 1919, Jeffcott [2] identified 
the critical speed, not as a speed of instability, but 
as a natural frequency of lateral vibration. By includ- 
ing damping in his analysis, Jeffcott showed that stable 
operation was possible above the critical speed. 

The increased demands on rotor systems for more 
power and lighter weight soon led manufacturers to 
design rotors for operation above the first critical 
speed. This led to instability problems, as rotors 
reached speeds above which stable operation was impos- 
sible. Several mechanisms of rotor instability have 
been identified and are the subject of works by 
Ehrich [3,4]. These include hydrodynamic bearings and 

seals, internal rotor friction, and turbine blade 
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effects. Each of the mechanisms of instability possess 
a common characteristic in that they produce forces 
tangential to radial deflections. This cross-coupl ing 
effect generally increases with speed until it overcomes 
the restoring and dissipative forces and the rotor 
becomes se if -exc ited . Violent whirling or "whipping” of 
the rotor occurs as the stability threshold is reached. 

While most instability mechanisms are nonlinear 
phenomena, they can generally be modelled as linear 
systems with variable coefficients. A significant non- 
linearity occurs, however, in bearings with radial 
clearances or "deadbands". Although these are not 
usually considered mechanisms of instability, their 
effect is significant on rotor performance. This effect 
has been examined by Childs [5] and others [6,7] whose 
results show that these nonlinearities often result in 
subsynchronous oscillations at frequencies exactly 
one-half of rotor spin speed. It has also been shown 
that bearings with radial clearances do not affect the 
overall stability of rotor systems. 

The need for still greater performance from rotor 
systems requires operation at speeds above what is now 
the limit of stability. For this reason, the focus of 
much recent work has been on the improvement of rotor 
stability characterist ics . This area has appeared more 
promising with the advent of active magnetic bearings. 
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Hab e pro a n d[iu L i aiu [ 8 ] P p e s en t an sxcs 1 1 6nt discussion 
on the practical applications of these bearings. 
Schweitzer and Lange [9] present a more detailed 
description of magnetic bearing characteristics, and 
show how they can be implemented in a closed-loop 
control system. GondhaleRar and Holmes [10] discuss the 
various designs of magnetic bearings for controlling 
rotors. They show, both analytically and experimen- 
tally, how suitable configurations of magnetic bearing 
systems can lead to linear, uncoupled, completely 
determined control forces. 

For stabilizing rotor systems, Burrows and 
Sahinkaya [11] examine the use of magnetic bearings to 
control oil-whirl. They apply pole-placement techniques 
to a single-mass, rigid rotor supported by hydrodynamic 
journal bearings. They determine the contol gains 
neccessary to position the unstable pair of eigenvalues 
sufficiently far from the imaginary axis, and they 
discuss the constraints which affect reachable pole 
zones. More recently, Stanway and Burrows [12] examine 
rotor stabilization by applying control to the rotor’s 
support structure. They find that, with full state 
feedback, it is possible to obtain some degree of stabi- 
lization. They discuss the use of observers to construct 
the full state vector when it is not directly obtain- 
able. Schweitzer [13] examines stabilization, using 
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magnetic bearings, of a linear multi-body rotor with 
internal damping. He uses one magnetic bearing set and 
defines an approach to locate the bearing based on the 
rotor’s mode shapes. He determines an "optimal” control 
by maximizing the lowest damping of the system’s modes. 
In a later work [14] he applies magnetic bearing control 
to a low order model obtained from a higher order modal- 
based model. Schweitzer then develops the control for 
the lower order model and examines the "spillover" 
effects on the higher order system. 

The purpose of this study is to investigate the 
stabilization of self-excited rotor systems by active 
control. In particular, the application of linear 
control to nonlinear rotor systems under the influence 
of bearing deadband is examined. Its importance results 
from the fact that nearly all physical rotors supported 
by radial rolling element bearings are affected by 
deadband. While this factor is usually neglected in 
most rotor studies, it is one which must be understood 
for the full benefits of stability enhancement to be 


achieved . 



CHAPTER III 


THE TWO DEGREE -OF -FREEDOM MODEL 

Introduction 

The first phase of this study deals with the 
stabilization of a two degree-of -f reedom model which is 
similar to a modified Jeffcott model. The equations of 
motion for the model are developed and converted to 
dimensionless form. The control strategy for rotor 
stabilization is developed for the linear rotor and 
applied to the linear and nonlinear models. Maximum 
control bearing forces required to stabilize the rotor 
are taken from dynamic simulations of linear and nonlin- 
ear models. Finally, the frequency spectrum of the 
nonlinear simulation is analyzed to observe the charac- 
teristic behavior of the nonlinear system. 
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Fig. 3.1. Two Degree-of-Freedom Model. 
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of the types of behavior found in more complex rotor 
systems . 

As shown in Fig. 3.1, the model consists of a uni- 
form disk of mass m suspended by a massless rigid shaft. 
The bearings have a radial clearance g which provides a 
"deadband" in the force-deflection curve. The shaft is 
supported on each end by massless radial bearings having 
a combined linear stiffness K&. Damping is introduced 
by hydrodynamic seal type elements which provide stiff- 
ness and damping K s and C s , respectively, and cross- 
coupled stiffness and damping Q s and Cq, respectively. 
Although for true rotor systems most parameter values 
are functions of rotor spin speed, the most significant 
is the cross-coupled stiffness Q s . Since this term is a 
source of self-excited vibration, it is sufficient to 
introduce speed dependency with this term only. Thus, 
it will be given the familiar approximate form for 
hydrodynamic seals, 

C s w 

Q s = . (3.1) 

2 

The rotor is unbalanced with the mass m at a radius e. 
The rotor is assumed to be symmetric, both radially and 
laterally. The rotor spin speed u> is constant, although 
different speeds will be examined. Motion occurs in the 
Y and Z directions only. 
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A force diagram is shown in Fig. 3.2, where Ft, z 
and Ffcy are the restoring forces from the support bear- 
ings, F s2 and F S y are the "seal" forces, and mew 2 is the 
unbalance force. Fig. 3.3 shows a force-deflection 
curve for the support bearings. The magnitude of the 
radial force produced by the bearings is given by 


F br = 


0 , | r | < g 

K (r - g) , Irl > g 
b 1 1 


(3.2) 


with the force acting in the direction of radial 
displacement. Resolving this force into the Y and Z 
directions yields 


F by = 


, | r j = J y* + 2 2 < g 


(3.3 ) 


K ( i - — )y , I r I > g 

» | r | 


and 


F b2 


0 , | r | < g 

6 

K (1 - - )z , | r J > g 

b I r» I 


(3.4) 


The forces due to the seals are given by 

F sy = K s y + C s y + Q s z + C q 2 (3.5) 


and 


F s 2 : K s z + c s^ " Oi S y - C q y . 


(3.6) 



Support Bearing Force-deflection Curve 
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Combining the force relationships in their respec- 
tive directions yields the equations of motion for the 
two degree-of -freedom rotor. 

8 

my + C s y + K b (l )y + K s y + C q z + Q s z 

| r I 

= meu> 2 cos (wt ) (3.7) 

« 

mz + C s z + K b (l - - )z + K s z - C q y - Q s y 

I r I 

= meu> 2 s in (wt ) (3.8) 

Division by m yields 

.. C s K b 8 K s C q __ Q s 

y + — y + — (1 - — )y + — y + — z + — z 

m m | r | m m m 

= eu> 2 cos (wt ) (3.9) 

and 

.. C s . K b g Kg C q ^ Q s 

2 + — z + — (1 - — )z + — z - — y - — y 

m m | r | m m m 

: ew 2 s in (wt ) . (3.10) 

To effect a more general study, these equations 
are nond imens ional ized . First, a system frequency is 
defined as the undamped, uncoupled frequency given by 

K s + K b 
w o d = 


m 


(3.11) 
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By defining the dimensionless time parameter r as 

T = w 0 t (3.12) 

and using the chain rule it follows that 


dy 

dy 

dr 

dy 


y = : 

: 

: 

: w 0 — = w 0 y' . 

(3.13) 

dt 

dr 

dt 

dt 



Simi larly, 


y = 5 0 2 y". 

2 = W 0 Z # , 


and 


(3. 14) 

(3.15) 


2 = So 2 *"- (3.16) 

Substituting for t, using (3.13) - (3.16) and dividing 
by w 0 2 yields 



and 
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Z ** + Z * + 

mw 


*b 8 K s 

(1 - - )z + z 

mw 2 i r . mw 2 

0 11 o 


Y* 


mw 


mw 0 2 


w c 

y = — s in { — 



<: 


Now make the following definitions: 


w 

- = R, 

Wft 


Ss = 


2mw, 


Sq = 


2mw, 


R kb = 


K b 


K g + K b mw 0 2 


Rrs = 


K s + K b mw 0 2 


C s w Q s c s w 

Q s = • so = = S s R . 

2 mw 0 2 2mw 0 w 0 


Substituting these yields 


y" + 2C s y' + R kb ( 1 - - )y + Rk S y + 2 Sqz 

I r I 


+ ?gRZ = R 2 ecos(Rr) (3 


and 


z" + 2? s z' + R kb (i - - )z + R ks z - 2C q y' 


- C s Ry = R 2 esin(Rr). (3 


Finally, division by a characteristic length, 


. 18 ) 


• 19 ) 

. 20 ) 
say 


e, yields the d imens ion 1 es s model 
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y " + 25 s y # + RKbd 


and 


8 

- - )y + Rfc S y + 2£ q Z' 


+ £ s Rz = R 2 cos(Rt) 


(3.21 ) 


8 

z" + 2S S Z' + RRbf 1 - - )Z + Rr S 2 - 25q7' 


- 5 s Ry = R 2 s in (Rr ) . (3.22) 


Table 3.1 summarizes the dimensionless parameter defini- 
tions . 

Equations (3.21) and (3.22) represent the non- 
linear, uncontrolled rotor. Defining the state vector x 
as 


x = [y z v y v z ] T 

the system in first-order form becomes 
x * z Ax + d 

where 


0 

0 

l 

0 

0 

0 

0 

i 

(R + R ) 

-c R 

-25 

-25 

Ks kb 

s 

s 


5 R 

-(R + R ) 

25 

-25 

s 

ks kb 

q 



(3.23) 


(3 . 24) 


(3.25) 
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and 


d 


0 

0 


R 2 cos (Rt ) + 
R 2 sin(Rt) + 


R 

Kb 


g 

- )Y 


R 

Kb 


g 

- )z 


(3.26) 


Note that the vector d contains both the unbalance 
forces and the nonlinear contributions from the support 
bearings. The same restrictions apply as in equations 
(3.3) and (3.4) to the nonlinear terms, i.e., 

I 

R Kb(l - = ) : 0, If I F 1 < g . 

I r l ' ' 

Also note that, by definition, Rj^s + Rub = l, so the 
system matrix becomes 


0 

0 

1 

0 

0 

0 

0 

i 

-1 

-c R 
s 

CO 

OJ 

1 

-25 

q 

5 R 
s 

-1 

25 

q 

-25 

s 


Dimens ionless 


i Nomenclature 

Frequency 

Damping 
St if f ness 

Cross-Coupled 
Damping 

Cross-Coupled 
St if f ness 

Spin Speed 

Displacements 

Deadband 


Time 
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trameter 


Def init ions 
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Characteristics of the Uncontrol 1 ed Rotor 

This study concerns stabilization of rotors for 
speeds at which the uncontrolled rotor is unstable. For 
nonlinear systems, however, stability and instability 
are amplitude dependent. When examining the instability 
of a rotor with deadband, it is apparent that the effect 
of the deadband is diminished as the rotor’s amplitudes 
of vibration begin to grow. Global stability can be 
determined by neglecting the deadband, producing a lin- 
ear system compatible with various linear analysis tech- 
niques. Therefore, this study deals with analyzing lin- 
ear systems obtained by neglecting deadband. The 
results of these analyses are applied to the nonlinear 
system, which is examined by simulation. 

The eigenvalues of the linear system can be 
readily obtained and will provide nearly all of the 
information required, both for examining the stability 
of the rotor and for determining the desired control 
laws. Examination of the system matrix A shows a speed 
dependency in the term £ S R. Thus, for any given speed 
R, a different linear system is obtained. At some 
speed, the system becomes unstable. That speed is the 
onset speed of instability or the instability threshold 
speed. Using the parameter values shown in Table 3.2, 
the system eigenvalues over the speed range 0 < R < 5 
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are shown in Fig. 3.4. The system is at the threshold 
of instability for a speed of R = 2.2. The frequency of 
the unstable mode ("whirl") is 1.1, approximately one- 
half that of the rotor spin speed, characteristic of 
self-excited rotors. As expected, the nonlinear system 
with deadband becomes unstable above R = 2.2. This is 
shown by simulation results depicted in Fig. 3.5, using 
a deadband g = 1 for R = 2.3. 


Table 3.2 

System Parameter Values 


Parameter 


Kb 


Value 
0.5 
0 . 1 
0.8 
0.2 


UNCONTROLLED ROTOR POLES 



Uncontrolled Rotor Pole Trace. 
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ROTOR CENTER POSITION 


Uncontrolled Rotor at R = 2.3 



"2.5 -1.5 -0.5 0.5 1.5 2.5 

(Thousonds) 

Y 


Fig. 3.5. Unstable Response of Nonlinear System at R * 2.3. 
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The Contro 1 Actuator 

The problems of applying completely definable 
forces to high speed rotating machinery have been the 
focus of much recent study [8, 9, i 0 ] . The result of this 
effort is a versatile class of devices widely known as 
magnetic bearings. These bearings utilize electromag- 
netic interaction to produce forces on shafts suspended 
within their housings. Many experimental studies and 
several practical applications attest to the feasibility 
of these bearings for use in rotating machinery. 

The force produced by a magnetic bearing can be 
made relatively proportional to the driving current. The 
use of two bearings grouped together results in a device 
which can produce two independent forces in arbitrary 
directions. By structuring the input currents to be 
functions of the states of the system, a closed-loop 
control system is produced. The forces can be made to 
resemble stiffness, damping, combinations of these, or 
other suitable forms. 

For the purposes of this study, the magnetic bear- 
ing forces are assumed to act in orthogonal directions 
coinciding with the Y and Z rotor coordinate directions. 
The transients associated with the actuator and control 
loop are assumed to be negligible compared with those of 
the rotor. The input current is proportioned to the 
state vector x through a gain matrix G so that 
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u 


u 

1 
u 

2 


= -Gx 


where 


(3.28) 



g, 

g . „ 

8 

g 

G = 

1 1 

12 

13 

14 


g„ 

g 

8 

g 


21 

22 

23 

24 


(3.29) 


Stabil izing the Rotor 

Modifying equation (3.24) to include the control 
forces acting through a control distribution matrix B 
yields 


X * = Ax + Bu + d . (3.30) 

Again, since system stability is of interest, 
bearing deadband is temporarily neglected. The stability 
of linear systems is independent of input, so the vec- 
tor d is neglected temporarily as well. The resulting 
1 inear system is 


X ' : AX + B U 


(3.31 ) 


or 


(A - BG)x. 


(3.32) 
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Since there are independent, unconstrained control 
forces acting in each direction, it can be shown that 
all states are accessible by the control, and the system 
is controllable. Thus, by selection of the gain matrix 
G, the eigenvalues of (A - BG ) can be arbitrarily posi- 
tioned . 

A simple method for stabili2ing the rotor can be 
developed by examination of the closed-loop system mat- 
rix (A - BG ) . From the definition of the state vector x 
and the form of the input u, the control distribution 
matrix B has the form 



(3.33) 


This results in the closed-loop system matrix 


(A - BG) = 


0 

0 


0 

0 


- ( 1 + g ) ~(5 R+g ) 

11 S 12 

(5 R~g ) -(l+g ) 

S 21 22 


1 

0 


-<2S + g 4 J 

S 13 

(25 -g ) 
q 23 


0 

1 

- (25 +g ) 
q 14 

-(25 +g ) 
S 24 


(3.34) 


Since the rotor is radially symmetric, it follows 


that the control should conform to the same symmetry. 
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Theref ore 

§21 = -gl2 

§22 = g21 

§23 = -gl4 

and 


(3.35) 

(3.36) 

(3.37) 


§24 = g 1 3 . 


Thus, 


(3.38) 


(A - BG) 


0 0 
0 0 


-(i+g ) -(? R+g 

11 S 12 

(5 R+g ) -(!+§ ) 

. S 12 11 


1 0 
0 1 


-(25 +g 

. ) 

"(25 + g 

) 

s 

13 

q 

14 

(25 +g 

„ ) 

-(25 + g 

„ ) 

q 

14 

s 

13 J 


(3.39) 

Examination of this matrix shows a speed depen- 
dency only with the term (5 S R + gig). By selection of 
gig as a linear function of R, this term can be held 
constant for any spin speed. Furthermore, by setting 
§11 : §13 = §14 = °i the above matrix at an arbitrary 
spin speed can be made equal to the open-loop matrix A 
at some speed, say R 0 . In other words, making 


§12 ~ Cs(^o ~ ® ) 


(3.40) 


will make 
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(A - BG) R = A Ro . (3.41) 

This essentially cancels the effect of the cross-coupled 
stiffness, which is the mechanism of instability in this 
model. Furthermore, the eigenvalues of (A - BG ) at any 
spin speed R remain constant and equal to those of A at 
the speed R 0 . Thus, the problem of control determina- 
tion is reduced to choosing the speed R 0 for which the 
system has the desired characteristics. 

One logical choice of the "reference speed" R 0 can 
be made by examining the pole-trace of Fig. 3.4. At the 
speed R = 0, the system has the highest relative damping 
of any speed. Relative damping associated with the i th 
mode is defined as 

-«i 

tRD * / ; (3.42) 

J <5i* + u> ^ 

where X* = + Jwj represents the i th eigenvalue of the 

system. As an eigenvalue moves towards the unstable 
zone, its relative damping decreases. The minimum rela- 
tive damping thus becomes a useful measure of the sta- 
bility of the system. At R s 0, this value is maximum 
(Crd = 0.497), with both modes being equally damped. 

For comparison, the minimum relative damping at the 
uncontrolled first critical speed of R = 0.78 is 
S R D = 0.323. Thus, the choice of R 0 = 


0 results in an 
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"optimally” stable system with a 54* increase in criti- 
cal speed damping. ("Optimal" of course, refers only to 
the control approach chosen here.) 

Using equation (3.40) with R 0 = 0 yields 

6 1 2 = -S s R (3.43) 

and 


0 -5 R 0 

s 

5 R 0 0 

s 


0 

0 


(3.44) 


This Keeps the poles of the closed-loop system at 
Xj = -0.55723 + J(0.97366), and 
X 2 : -0.44272 + j(0.77366) 
for any rotor spin speed R. 

Contro 1 Performance 

While knowledge of the system poles gives an indi- 
cation of the degree of stability and damping of the 
system, it gives no insight into other factors relative 
to the performance of the control. For example, while 
the implemented gains may cause a stable eigenvalue 
assignment, the force levels required to accomplish the 
stabilization may be higher than the physical actuators 
can produce. Also, the control of the nonlinear system 
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is of concern. These questions can be addressed by 
dynamically simulating the controlled rotor. 

Using the control developed in the previous sec- 
tion, numerous simulations were conducted for various 
spin speeds and bearing deadband values. Fig. 3.6 shows 
a typical simulation result. For each case, the maximum 
control force required during the steady-state rotor 
response was computed. These are plotted in Fig. 3.7 
versus spin speed. The three cases include the linear 
system (g = 0) and two nonlinear systems (g = 1, and 
g = 10). This figure shows a nearly linear dependence 
of the maximum control force on the spin speed R. 
Furthermore, the variation between linear and nonlinear 
cases is slight, if not negligible. Another interesting 
performance indicator is how well the rotor responds to 
an impulsive disturbance. Fig. 3.8 through Fig. 3.10 
show the rotor response to an impulse velocity in one 
direction. Again, the cases include the linear system 
and two nonlinear systems, showing the maximum control 
force for each. While for these cases the maximum 
forces are higher and increase with deadband, this is to 
be expected. Loss of stiffness due to deadband causes a 
greater excursion from the steady-state orbit, producing 
a higher control force. 

Finally, the frequencies of a rotor orbit are 
often of interest. Sub-synchronous whirls cause cyclic 
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flexure of the rotor shaft and are often destructive. 
They also indicate the presence of self-excited vibra- 
tion mechanisms which lead to instability as the speed 
is increased. Fig. 3.11 shows the frequency composition 
of the rotor response of Fig. 3.10. it is predominantly 
synchronous, as the se 1 f -exc it ing forces are counter- 
acted by the control, and subsynchronous components are 
not detectable. 

Summary 

In this chapter, a two degree-of freedom model was 
developed and stabilized for an arbitrary speed . While 
the control approach utilized was fairly simplistic, it 
was shown to be quite effective. More importantly, it 
allowed for an evaluation of the control performance. 
While such indicators as control force levels and the 
response of the nonlinear rotor may change quantita- 
tively for different control schemes, the trend of each 
would probably not differ substantially. The results of 
this chapter will be applied to a more realistic, com- 
plex rotor system in chapter four. 
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Fig. 3.7. Maximum Required Control Force Vs. Spin Speed. 
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Fig. 3.10. Rotor Response to Impulse at R * 5, g = 10. 





CHAPTER IV 


THE SIX DEGREE-OF-FREEDOM MODEL 


Introduction 

This chapter deals with stabilizing a six degree- 
of-freedom rotor using one active magnetic bearing set. 

In the previous chapter, a two degree-of -f reedom rotor 
was stabilized by using the control bearing to counteract 
se If -exc it ing forces. The result was a rotor that was 
stable for any given spin speed. In many physical 
rotors, however, the locations of se If -exc it ing mecha- 
nisms dimensionally restrict the positioning and number 
of control forces that may be applied. Full state feed- 
back for large rotor systems would also be impractical. 
For these reasons, the six degree-of -freedom rotor model 
is examined. This model enables the introduction of 
several self -exc iting mechanisms and nonlinearities 
located away from the point of control application. To 
investigate the stabilization of this rotor, the model 
is presented in d imens ion 1 ess form. Its uncontrolled 
characteristics are studied to determine such factors as 
control location and signal measurement. An "optimal” 
stability criterion is defined using the eigenvalues of 
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the linear system. The control gains are determined by 
a search method. Finally, simulations are conducted to 
determine required control force levels for the linear 
and nonlinear systems. 

The Mode 1 

The six degree-of -f reedom model is shown in 
Fig. 4.1. This rotor has three uniform lumped masses 
m*, m 2 , and m 3 connected by two massless, flexible 

shafts of stiffness Kfi and Kf£* Acting at each mass 
are support bearings Kj, K 2 and K 3 , respectively. These 
are radial type bearings, each with a deadband g^, 

(i r 1 , 2, 3). Also at each mass are elements which pro- 

vide cross-coupled stiffness and damping Q s * and c qi- 
as well as support stiffness and damping K s ± and C si . 
Each mass is "unbalanced" by its mass m^ at a radius 
where all three unbalanced radii are assumed to lie in 
the same plane. The rotor is assumed to be radially 
symmetric operating at a constant spin speed w. Only 
motion in the Y and 1 directions is considered. Gyro- 
scopic effects, which can be shown to enhance rotor sta- 
bility, are neglected. 

The support bearings are initially assumed to have 
zero deadband. This simplifies the form of the equa- 
tions of motion. The nonlinear contributions from dead- 
band are included during the simulations by a 
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Fig. 4.1. Six Degree-of-Freedom Model. 
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disturbance vector. The resulting system is linear with 


the following equations of motion: 

m l 7 i + Cjyi + (K s i + Kjjyj + K fl (y t - y 2 ) 

+ CqiZj + QsiZj = miw 2 e 1 cos (wt ) (4.1) 

miz’i + CiZ! + (K s i + Kilzj + KfiUi - z 2 ) 

- - Q s iyi : m 1 u> 2 e 1 sin(wt ) (4.2) 

m 2 y 2 + c 2 y 2 + (Kg 2 + K 2 )y 2 + K fl (y 2 - yi ) 

+ K f2 (y 2 - y 3 ) + C q2 z 2 + Q s2 z 2 = m 2 w 2 e 2 cos (wt ) (4.3) 

m 2 z 2 + C 2 z 2 + (K s2 + K 2 )z 2 + Kf t (z 2 - Zj ) 

+ K f2 (z 2 - z 3 ) - C q2 y 2 - Q s2 y 2 = m 2 u> 2 e 2 s in (u>t ) (4.4) 

m 3 y 3 + C 3 y 3 + (K s3 + K 3 )y 3 + K f2 (y 3 - y 2 ) 

+ C q3 z 3 + QS3z 3 = m 3 u>2 €3COS (wt } (4.5) 

m 3 2 *3 + C 3 Z 3 + < K s3 + k 3) z 3 + K f2< z 3 ’ z 2 ) 

' Cq3Y3 ' Q S 3Y3 = rn 3 w 2 e 3 sin(wt ) (4.6) 


By defining the dimensionless parameters shown in 
Table 4.1, the equations are nondimens ional ized and 
become: 


R kl B fl 

yi" + 2SiYi' + 7i + (yi - y 2 ) 

R ml R ml 


R ml 


(4.7 ) 
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R K 1 R fi 

z l" + 2ClZl' + Zl + ( 2 i - Z 2 ) 

R ml 

Sl R 

" 2 5qiyi # “ ?1 = R el R 2 sin(Rr) (4.8) 

R ml 


_ r K 2 R f 1 Rf2 

y 2 " + 25 2 y 2 ' + y 2 + (y 2 - y t ) + (y 2 - y 3 ) 

R m2 R m2 Rm2 

C 2 r 

+ 2 £q2* 2 ' + 2 2 = R e 2 R 2 cos (R t) (4.9) 

R m 2 


r K 2 R fi Rf 2 

z 2 " + 2S 2 Z 2 ' + z 2 + (Z 2 - Z t ) + (z 2 - z 3 ) 

R m 2 R m 2 Rm 2 

5 2 r 

* 2 Cq 2 y 2 # - y 2 = R e 2 R 2 sin(Rr) (4.10) 

R m 2 

r K 3 R f 2 

y3" + 2 C 3 y 3 ' + y 3 + ( y 3 - y 2 ) 

R m3 R m3 


53 r 


+ 2 ^q3 Z 3 # 

- w 

23 

= R e 3 R 2 cos(Rr) 

(4.11) 

R m3 




r K3 

R f 2 




Z 3 " + 2 C 3 z 3 ' + z 3 + 

(z 3 - 

z 2 > 


R m3 

R m3 



' 2 ?q3y3 # 

5 3 r 


= R e 3 R2sin(Rr) 


- 

Y3 

(4.12) 


R ml 





By choosing the state vector x to be 


2 ? = [7l ?2 ?3 v yl v y2 v y3 z t z 2 z 3 v zl v z 2 v z3 )T (4.13) 
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the system in first-order form becomes 


x # = Ax + d 


( 4 . 14 ) 


where 


(R + R 


+ R + R 


( 4 . 15 ) 


( 4 . 16 ) 


( 4 . 17 ) 


0 


0 
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R 

ml r R 



R 

m3 


(4. 18 ) 


H 


2 £ 0 0 

qi 

0 25 0 

q2 

0 0 25 

q3 J 


(4.19) 


and 

d = disturbance vector of unbalance and nonlinear 
bearing forces. 


Characteristics of the Uncontro 1 1 ed Rotor 

This section describes the characteristics of the 
rotor model developed in the previous section. Table 
4.2 shows the parameter values used for this model. 

As in the previous chapter, stability can be 
determined by examination of the poles of the linear 
system. Fig. 4.2 shows a trace of the rotor poles for 
0 < R <5. The first mode is unstable above R = 1.78, 
the second above R = 4.1 and the third above R = 4.65. 
Fig. 4.3, Fig. 4.4, and Fig. 4.5 show modal "orbit" 
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Table 4 . 1 


Dimensionless Parameter Definitions 




Cross-Coupled 

Damping 


Cross-Coup 1 ed 
St if f ness 


Spin Speed 
Displacements 

Deadband 




Where: = m^ + m2 + m3, £t = e l + e 2 + e 3> and 


Kt = Kj + K s j + Kg + Kgg + K3 + K53 
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shapes for each of these modes at their respective 
unstable speeds. For the chosen mode, the modal orbits 
are obtained by combining the complex eigenvector and 
eigenvalue into the form of a solution and retaining the 
real part. These orbits show that the (mainly) rigid 
body “bounce" mode is the first to become unstable, fol- 
lowed by the rigid body "pitch” mode, and finally the 
flexible mode. This is the expected order, as the lower 
frequency (rigid) modes are most easily excited. 
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Table 4.2 

Rotor Parameter Values 


i 

i 

i 

i 


Parameter 

Value 

R 

0 . 25 

ml 


R 

0 . 50 

m2 


R 

0.25 

m3 


5 . 

0.25 

5 

0 

2 


5 3 

0 . 25 

R 

0 . 50 

K 1 


R 

0 

K2 


R 

0 . 50 

K3 


R 

0 . 50 

f 1 


R 

0 . 50 

f 2 


c 

0 . 05 

qi 


i 

0 

q2 


l 

0.05 

q3 


R 

0 

ei 


R 

1 . 00 

e2 


R 

0 

£3 




UNCONTROLLED ROTOR POLES 
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Fig. 4.2. Uncontrolled Rotor Pole Trace. 



MODAL ORBIT SHAPE 

First Unstable Mode at R — 1.775 
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Fig. 4.3. Modal Orbit Shape of First Unstable Mode at R = 1.775. 



MODAL ORBIT SHAPE 

Second Unstable Mode at R — 4.101 
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Fig. 4.4. Modal Orbit Shape of Second Unstable Mode at R = 4.101. 



MODAL ORBIT SHAPE 

Third Unstable Mode at R -» 4.65 
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Fig. 4.5. Modal Orbit Shape of Third Unstable Mode at R = 4.65. 
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Contro 1 Structure 

In this section, the control structure is devel- 
oped and implemented on the linear rotor system. 

I 

I 

Disturbances are neglected. Taking the system without 
disturbances, the control u is added so that 

x ' = Ax + Bu . (4.20) 

Now let 

u = -Gy (4.21 ) 

where 

y = Cx. (4.22) 


The above control form is chosen because full state 
feedback is impractical on physical rotor systems. 

The closed-loop system becomes 

X' = (A - BGC )x . (4.23) 

To produce the two control forces 

(4.24) 

the gain matrix must be based on the form of the 
output y. The output y is chosen to take the form 



where 


1 = C Y c V yc Z c V zc ]T 


5 i 

(4.25) 


Y c is a linear combination of the Y-direction 
displacements y^, y 2 , y3, 

Vy C is a linear combination of the Y-direction 
velocities v y i, v y2 , v y3 , 

Z c is a linear combination of the Z-direction 
displacements 2j, 2 2 , 23 , and 

V zc is a linear combination of the Z-direction 
velocities v zl , v z2 , v z3 . 


Obviously, there are many ways to structure the output 
2 - The above choice is logical, however, in that it 
simplifies the structure of both the output matrix C, 
and the resulting gain matrix G. Recalling the radial 
symmetry of the rotor, it is apparent that the Y and Z 
directional measurements of the output should be symmet- 
ric. This results in C having the form 



0 


b 

0 

0 

0 


c 

0 

0 

0 


0 

d 

0 

0 


0 0 0 0 0 

e f 0 0 0 

0 0 a b c 

0 0 0 0 0 


0 0 
0 0 
0 0 
d e 



(4. 26) 
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Also from symmetry and the results of the previous 
chapter, the gain matrix G has the form 



g 

g 

g 

g 

“ 

il 

12 

13 

14 


-g 

-g 

g 

g „ 


L 13 

14 

1 1 

12 J 


The apparent difference in equation (4.27) and the gain 
matrix of the previous chapter results from the order of 
the state vector. The control distribution matrix B 
depends upon the chosen point of application of the con- 
trol forces. It is assumed that both forces will be 
applied at one mass only. 

Contro 1 Determination 

To stabilize the rotor represented by the closed- 
loop system of equation (4.23), it is sufficient to 
position all of the poles in the left half plane. It is 
not necessary, however, to place all the poles to 
achieve stability. Thus, the system of equation (4.23) 
need not be completely controllable. By using one mag- 
netic bearing set in the absence of full state feedback, 
the control is suboptimal. 

The various analytical techniques available for 
dealing with the linear suboptimal control problem gen- 
erally influence a portion of the system in a specific 
way, and minimize the influence of the control on the 
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remainder of the system. Modal control, for example, 
"fine-tunes" a specific mode or modes, while the "spill- 
over" of the control to other modes is not directly 
controlled. The use of this technique for this applica- 
tion limits the number of directly controllable modes to 
one. While this allows the rotor to operate above 
previously unstable speeds, the influence of the control 
on the two remaining modes would be questionable. 

For these reasons, a direct search method for 
determining the control gains is used. Using Powell’s 
method [15] a function is minimized by iteratively 
searching an independent vector, in this case, the con- 
trol gains. As discussed in the previous chapter, a 
useful measure of the stability of the system is the 
relative damping of the systems’s modes. The relative 
damping Sup of the i*- 1 " 1 mode is defined as before to be 


Srd 1 



(4.28) 


where Ai : <5j + is the eigenvalue associated with 

that mode. It is the minimum relative damping which is 
critical to the stability of the system. Maximizing 
that minimum relative damping is equivalent to maxi- 
mizing the shaded region of Fig. 4.6. The minimum rela- 
tive damping of the closed-loop system is determined by 
a computational function. The negative of this function 
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is minimized using Powell's method to determine the con 
trol gains. The applications oi this approach are 41s- 

cussed in the next section. 
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Contro 1 AppI icat ion 

Before the control gains are determined, two 
factors must be specified concerning the control: 

(1) Point of application of the bearing, 

(2) Combination of measurements for the output 
s ignal . 

To best determine the location for applying the 


control, it is useful to examine the uncontrolled 
rotor’s modal orbit shapes. Recalling Fig. 4.3 through 
Fig. 4.5, it can be seen that, for affecting the first 
unstable mode (Fig. 4.3), the bearing can be located at 
either three mass positions. This is because that mode 
is very nearly a rigid body translational (or "bounce”) 
mode. Examination of the next unstable mode (Fig. 4.4) 


clearly shows that a bearing located at the center mass 
location offers no stabilizing effect for this mode. 
Since this mode is unstable at a speed lower than the 
third mode, controlling it takes priority over modes 
that follow. Therefore, the magnetic bearing should be 
located at an end mass location, say that of 114 . 

The decision of where to measure and how to com- 


bine the signals for the output % is not as simple as 
the location of the control. It is helpful to make a 


restriction regarding this factor. Since full state 
feedback is impractical for most rotor systems, it is 
assumed that the information from all three mass 
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locations will not be available. The output can contain 
information from two locations, at most. An examination 
of the uncontrolled rotor’s modal orbit shapes shows 
that if information is provided from the center mass 
location only, the second ("pitch") mode cannot be 
affected. Thus, an end location must be measured. For 
stabilization of the third mode, information from the 
center mass is required as well. To gain more insight 
into the question of measurement, the search program was 
run for varying combinations of measuring locations. 

For each trial, the stabilization capabilities were 
observed. The results of these show that the best degree 
of stabilization is obtained from signal measurement at 
the center mass location and the end location where the 
magnetic bearing is to be located. Furthermore, by 
appropriate weighting of the measurements an even higher 
degree of stabiliztion is achieved. 

Using the end mass m^ as a location for the mag- 
netic bearing and measuring the displacement and veloc- 
ity signals at the end and center masses m^ and mg with 
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the search program was used to determine the control 
gains. For various spin speeds the gains were computed 
to maximize the relative damping. For the chosen output 
and bearing location, the rotor was stabilized for 
speeds up to approximately R = 7.2. Fig. 4.7 shows the 
control gains obtained over the speed range 0 < R < 8. 
This figure illustrates that for any given speed, there 
are numerous solutions which optimize the damping. 

Using the gains obtained from a prior speed as starting 
values and incrementing the speed slightly, it was pos- 
sible to "tracK" a particular solution over the speed 
range. The irregularities observable in the figure, 
however, show the difficulties associated with this 
approach. 

The poles of the controlled rotor system using the 
gains of Fig. 4.7 are shown in Fig. 4.8. Enlarged views 
of various sections are shown in Fig. 4.9 through 
Fig. 4.11. Due to the somewhat irregular nature of the 
control gain solutions, the behavior of the poles, espe- 
cially at lower speeds, is somewhat erratic. 

For comparison, the control gains from Fig. 4.7 at 
the stability threshold speed (R = 7.2) are used and 
held constant over the entire speed range. Fig. 4.12 
shows the rotor poles for these constant gains. As seen 
in the enlarged views of Fig. 4.13 through 4.15, the 
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poles are more well behaved than those obtained using 
speed dependent gains. More importantly, however, is 
the fact that the system is stable over the same speed 
range. A close comparison of Fig. 4.8 with Fig. 4.12 
reveals that the overall behavior of the poles does not 
differ significantly. This indicates that, at least for 
this model, there is little advantage to using speed 
dependent gains. Therefore, the rest of this chapter 
will deal with systems controlled by the constant gains 
obtained at R = 7.2 
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CONTROLLED ROTOR POLES 
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Fig. 4.8. Controlled Rotor Pole Trace with Speed Dependent Control. 
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CONTROLLED ROTOR POLES (0£R<:8) 

Enlarged View 
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Fig. 4.11. Controlled Rotor Pole Trace with Speed Dependent Gains, Section C Enlarged. 




CONTROLLED ROTOR POLES 

(0 * R £ 8) Constant Gains 



Fig. 4.12. Controlled Rotor Pole Trace with Constant Gains 









CONTROLLED ROTOR POLES 

(0*R£8). Constant Gains, Enlarged View 
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Fig. 4.14. Controlled Rotor Pole Trace with Constant Gains, Section B Enlarged. 





CONTROLLED ROTOR POLES 

(0<R£8), Constant Gains, Enlarged View 
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Fig. 4.15. Controlled Rotor Pole Trace with Constant Gains, Section C Enlarged. 





68 


Contro 1 Peri ormanc e 

The performance of the control is evaluated by 
examining the maximum control force required to achieve 
stabilization. This control force is obtained by simu- 
lating the rotor at steady state, or in a limit cycle if 
one exists. Both the linear and nonlinear rotors are 
simulated. A power spectral density of a nonlinear 
simulation reveals its frequency components. 

Fig. 4.16 shows the maximum control force required 
to stabilize the rotor over the speed range 0 < R < 6. 
Though the rotor is stable for speeds up to R = 7.2, the 
amplitudes of vibration become unrealistically large as 
this speed is approached. For the nonlinear systems, 
the deadband g is equal in both the end supporting bear- 
ings. Fig. 4.16 shows that, as deadband increases, the 
force required to stabilize the rotor does also. The 
slight peak in the force curves for g = 0 and g = 1 
occurs at a speed slightly less than R = l. This is 
approximately the same as the first critical speed. 
Comparison of Fig. 4.16 with the corresponding plot of 
the previous chapter (Fig. 3.7) shows that nearly the 
same level of force is required for this model with the 
exception of that for the g = 10 case, which for this 
model requires an order of magnitude greater force. 

A typical simulation of the nonlinear system with 
g = 1 is shown in Fig. 4.17 through Fig. 4.19. At R = 2 
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this system exhibits mainly synchronous motion. Increas- 
ing the speed to R = 6, as shown in Fig. 4.20 through 
Fig. 4.22 changes the system dramatically. While the 
center mass exhibits mainly synchronous motion, the end 
masses exhibit complex subsynchronous motions. This 
results from the self -exciting mechanisms being located 
at both end masses, and none at the center. In all 
cases, the amplitudes at the mass where the control is 
applied are significantly lower than at the other 
masses. Fig. 4.23 shows the frequency components of the 
end mass response indicating mainly subsynchronous 
motion. The response of the rotor with g : 10 at 
R = 6.0 can be seen on Fig. 4.24 through Fig. 4.26. 

High amplitude subsynchronous motion is clearly evident, 
resulting in high required control forces to maintain 
stabi 1 ity . 

Summary 

In this chapter, a six degree-of -freedom model was 
presented and stabilized to the maximum possible speed. 
The control was supplied by one magnetic bearing set, 
optimally located by examining the rotor’s modal orbit 
shapes. For the model parameters chosen, this location 
was found to be the end mass of the rotor. The output 
measurement combination for the highest degree of stabi- 
lization combined the signals from the center mass and 
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the end mass where the magnetic bearing was located. 

The control gains were determined by use of a search 
method which maximized the damping of the least damped 
mode. Speed dependent gains offered little improvement 
over constant gains for this model. All three natural 
modes could be stabilized for speeds up to R = 7.2, 
which resulted in nearly a 300Z increase in the thresh- 
old speed of the original rotor. The control forces 
required to achieve stabilization were determined for 
various spin speeds and bearing deadband values, and 
were found to increase with increasing speed and bearing 
deadband. Though the magnitudes of these forces grew 
relatively large as the threshold speed was approached, 
their magnitudes for the linear and low value (g = 1) 
deadband systems were comparable to those found for the 
two degr e e -of -f r eedom model for speeds up to three times 
that of the original threshold speed. 


MAXIMUM CONTROL FORCE Vs. SPIN SPEED 
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Fig. 4.16. Maximum Required Control Force Vs. Spin Speed. 
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ROTOR CENTER POSITION 
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Rotor Response at Mass 2 with R = 6, g 


Fig. 4.21. 
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ROTOR CENTER POSITION 


R = 6.0 , Moss 3 





POWER SPECTRAL DENSITY 

Response at Mass 1, R — 6.0, g — 1 
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Fig. 4.23. Power Spectral Density of Response of Fig. 4.20. 









80 


ROTOR CENTER POSITION 


R = 6.0 , Moss 2 



-8.00 -6.00 -4.00 -2.00 0.00 2.00 4.00 6.00 8.00 

Y 

Fig. 4.25. Rotor Response at Mass 2 with R = 6, g = 10. 
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ROTOR CENTER POSITION 


R = 6.0 , Moss 3 
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Fig. 4.26. Rotor Response at Mass 3 with R « 6, g = 10. 



CHAPTER V 


CONCLUSIONS 

Using active control, the stabilization of 
self-excited rotor systems has been investigated. Two 
different rotor models have been used, and linear 
and nonlinear cases of each examined. For each case, 
the applied control has enhanced the stability charac- 
teristics of the rotor significantly, allowing it to 
operate above otherwise unstable speeds. 

For the two degree-of -freedom model, the control 
was structured to counteract se If -exc iting forces. In 
the absence of these destabilizing forces, the rotor was 
stabilized for any spin speed. The control forces 
required to maintain stability were found to vary nearly 
linearly with spin speed, and were not significantly 
affected by bearing deadband. While this control 
approach seemed somewhat simplistic, it was shown to be 
quite effective. The versatile nature of magnetic bear- 
ings could allow implementation of such a control scheme 
in a physical rotor, especially where dimensions would 
permit locating a magnetic bearing near each major self- 
exciting mechanism. 
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The control approach used for the six degree-of- 
freedom model was shown to be effective at improving the 
stability of the rotor, although the improvement was 
limited. The system was stabilized for speeds up to 
nearly 300* of the original instability threshold speed. 
The significant result was the stabilization of each of 
the rotor’s three natural modes using one magnetic bear- 
ing. The placement of the magnetic bearing was such that 
it maximized the effect that the bearing could exert on 
each of the rotor’s natural modes. The feedback signal 
for the control was structured to be a composite output 
signal, with less than full state information. To offer 
the most system improvement, it was found that this sig- 
nal should contain information from at least the center 
mass location and end mass location where the bearing 
was located. The control forces required to stabilize 
the rotor were shown to increase nonlinearly with spin 
speed and deadband. However, the relative magnitudes of 
these forces were comparable to those found for the two 
degree -of -f reedom model. Obviously, the particular 
parameters chosen to describe the model affect the 
improvement capabilities of the control. However, in 
many cases even a small improvement in the stability 
threshold speed is worthwhile goal. 

Several interesting topics relating to this study 
remain to be addressed. First, the function used to 
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determine the optimal control gains for the six degree- 
of-freedom model in this study was somewhat poorly 
suited for iterative methods. The parameter being maxi 
mized, the minimum relative damping, was a value asso- 
ciated with only one eigenvalue. As the eigenvalues 
"moved” under the effect of the control, the one with 
which the relative damping was associated would switch. 
This led to there being many solutions for the optimal 
gains, and erratic behavior of the eigenvalues under 
speed dependent gain solutions. One possible research 
question would be to examine alternative optimization 
criteria, possibly based on independent eigenvalue 
excurs ions . 

Secondly, the accurate modeling of many of the 
physical phenomena affecting rotor systems is presently 
an issue of uncertainty. These often highly nonlinear 
effects are generally treated by variable coefficients, 
the values of which may vary significantly from one 
effort to another. Of interest then, would be the 
sensitivity of a control application to variations in 
the modeling parameters. 

Lastly, the dynamics associated with the control 
loop should be investigated as they affect high speed 
rotor systems. Delays from signal measurement, process 
ing, and control actuators could greatly affect the 
performance of the control on machinery of this type. 
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